Pyroelectric vector and related polarization effects 
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The principal fact that the pyroelectric vector P° typical of polar crystals is invisible in equilibrium 
has been verified. Nevertheless, the general formula for P° as a definite bulk quantity associated 
only with the inversion properties of a given charge distribution in the bulk has been obtained. Its 
application to ferroelectric perovskites BaTi03, PbTiOa, and KNb03 is discussed. 

PACS numbers: 77.70.+a, 61.50.Ah, 61.50.Ks, 77.80.Bh, 77.84.Dy, 02.10.Gd 



For many years the pyroelectric vector stays to be one 
of the most mysterious objects in solid state physics, 
though the phenomena of pyro- and ferroelectricity are 
classical 0, ■ The attempts to grasp it result in the 
concept that this value can be described within an adi- 
abatic evolution between characteristic states 0. The 
elegant realization of this approach for electron subsys- 
tems Q is based on the formalism of the Berry phase 
H in which case the integration over a closed contour in 
a parametric space arises and is reduced to a quantized 
flow through that contour. Many important results are 
obtained within this approach |fj, la, IjJ. However, an 
evolution route over equilibrium states [s], is a subtle 
problem characteristic of thermodynamic charge systems 
[4J, as discussed below. One more problem is that the 
pyroelectricity is a structural property of polar crystals 
0, 0| and so must be typical of every crystal struc- 
ture without inversion, notwithstanding is the descrip- 
tion classical or quantum. It implies that apart from the 
quantum Berry-phase approach 0, 0] , other more clas- 
sical treatments should exist. 

Here we propose a novel approach to this subject and 
represent the pyroelectric vector P° in the form of a 
general invariant specifying polar crystals. Expressed 
in terms of characteristic length parameters and effec- 
tive charges in the bulk, P° is defined as a function of 
state, without recourse to any evolution. As a result, it 
can be estimated with making use of classical charge dis- 
tributions, though a quantum-mechanical calculation of 
effective charges is eventually assumed. 

Interested in periodic charge distributions in the bulk, 
we first point out that any local polarization P as a bulk 
property associated with the crystal symmetry [H HolU^ ] 
and so independent of a small external field E should give 
the field contribution PE to the specific energy, but this 
contribution is absent in equilibrium. To visualize this 
fact, we note that the projection P z along any princi- 
pal axis z is determined by the effective one-dimensional 
charge distribution p(z) = p(z + c), with the period c ap- 
propriate to the case, which arises after integrating the 
initial charge distribution over normal cross-sections of a 
unit cell. The value of P z can then be written down as 



P, 



dg 



zp(z)dz + / (z + c)p(z)dz 
Jo 



0. (1) 



Here the fact that the uniform field component E z does 
not distinguish any fixed choice of the unit cell bound- 
aries in the bulk is taken into account. So, the expres- 
sion in square brackets describes the polarization of the 
unit-cell portion of p(z) shifted by g from any initial one, 
while the outer integral with the factor 1 /c accounts for 
the averaging over all shifts as equiprobable. After inter- 
changing the order of integration in the final result 
follows from charge neutrality of a unit cell. 

This issue follows the statement by Larmor 0, ^| that 
any bulk dipolar properties are not typical of crystals in 
equilibrium. Such states are associated with the thermo- 
dynamic limit with sample surfaces adjusted so 
as to exclude depolarization effects £|Uj|. Notwithstand- 
ing, if a finite sample is metastable because of unstable 
charged surface [IT), the surface charges always form a 
macroscopic dipole which interacts with an external field 
and, being proportional to a sample volume, can be for- 
mally represented as an additive composition of dipolcs 
in the bulk. Such dipoles have, without doubt, the bulk 
nature at the first moment of forming a steady-state sur- 
face 0, . But due to the disturbance of a bulk sym- 
metry near the surface, initial polarization effects do not 
recognize polar crystals |U3, lla] • To single out true bulk 
effects, the scheme of a sample placed between shorted 
condenser plates normal to the polar axis is usually ap- 
plied 0, Q . But the effect of P° is then also removed in 
accord with so that only increments in P° caused by 
changes in temperature or stresses are actually recorded 
via transient signals Q. 

Since P° is evanescent for observation on polar crys- 
tals in equilibrium 0, flfl llfll ] and the polarity reversal is 
impossible as an equilibrium physical process in pure py- 
roelectrics, the definition of P° as a bulk measure of po- 
larity associated with the inversion symmetry of a given 
state alone should be addressed. To gain insight into this 
subject, we first consider a simple Bravais lattice com- 
posed of point charges q. We assume that this lattice 
is immersed in a uniform neutralizing charge background 
and are interested in a dipolar moment p of the volume of 
a unit-cell parallelepiped specified by elementary trans- 
lations ai, a.2, and a3 and situated anyhow in the bulk. 
If the inner point charge is displaced by r from the par- 
allelepiped center initially, then it is easy to show that 



upon any displacement R of that parallelepiped through 
the bulk, with casting r — R in terms of {a^}, p(R) takes 
the form 

p(R) = qK , q(r - R) = g R (mod {qa k }), (2) 

provided that Ro belongs to the parallelepiped volume. 
This result is evident, for when a recurrent point charge 
leaves the unit-cell volume, the inner portion of back- 
ground begins to neutralize the next point charge enter- 
ing. So, p(R) is a periodic saw-tooth function of R with 
vertical jumps and with zero mean value in agreement 
with 0. Sites of this Bravais lattice and midpoints in 
between are its inversion points. Thus, a unit cell cen- 
tered on any of those inversion points is described by 
p = 0. In the cases of midpoints it implies that any 
of q on the unit-cell boundaries is to be evenly shared 
between adjacent unit cells HH3, in accord with 10). 

Every regular structure of point charges of different 
species may be cast in terms of simple Bravais sublat- 
tices introduced. To describe their polarization, the one- 
dimensional treatment is expedient along directions nor- 
mal to each couple of parallel faces of the foregoing ele- 
mentary parallelepiped so as to exclude multiple projec- 
tions of every particular Bravais sublattice on the axis 
under examination. These directions are defined by the 
unit vectors e k and by the periods d k = (a^e^), where 



ei 



a 2 x a 3 



[|a 2 | 2 |a 3 | 2 -(a 2 a 3 ) 2 ] 1 /2 



(3) 



and the transformation ei — > e 2 — » e 3 is cyclic. 

We begin with a combination of two Bravais sublat- 
tices, which is then specified by their projections per nor- 
mal section of an elementary parallelepiped on a given k- 
direction, as shown in Fig. QJa) , where the label k on all 
the parameters is omitted for brevity. The foregoing four 
symmetrical choices of the unit cell lead to four different 
values of the dipolar moment: 



Pi 



(I12 - hi) q2 



P2 



P3= il2?2, P4 = -£l2?l 

P3 = -hiq2, Pi = hiqi 



( I12 - hi)qi 
2 

at l\ 2 < ^21: 

at l\2 > hi- 



Also keeping in mind the following combination 
P1+P2 (I12 - hi)(qi - #2) 



(4) 
(5) 

(6) 



antisymmetric to inversion, we conclude that the general 
antisymmetric invariant containing all the motifs typical 
of QJk© takes the form 

Z(qi, 92; W2, hi) = Ii2hi(li2 - hi)qiq2(qi - 92)- (7) 

Each of six zeros of J7J excludes any polarization and 
this set of zeros is complete. The symmetry of {7J also 
implies that in the case at hand the inversion can be 
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FIG. 1: (a) Linear structure of two point-charge sublattices in 
a compensating uniform background shaded, with indicating 
the symmetric unit cell choices, (b) Unit cell '1' after re- 
moving the symmetric part of charges by means of inversion 
relative to its center, with q = 92/2. (c) The further modifi- 
cation of unit cell '1' allowing for the inversion antisymmetry 
relative to the unit-cell ends. 



defined in the four equivalent manners: (i) /i 2 — » — 1\% 
and hi -> ~hi\ (U) qi -qi and q 2 -> -q 2 ; (hi) 
I12 hi', (iv) qi <-> q 2 - Note that the transformation 
(<Zi; £12) *-> (92^21), that is the combination of (iii) and 
(iv), results in pi «-> p 2 , but p 3 — > P4 takes place after 
the additional shift T9i(^i2 + '21)1 where ?i2 + l 2 \ = 
in agreement with (J2J. Hence, those couples of cells are 
not mixed by operations related to inversion and so may 
be discussed separately. 

To create (JTJ normalized per unit cell, we start from 
unit cell T', where the antisymmetric charge combina- 
tion arises as a result of inversion relative to the unit-cell 
midpoint, as shown in Fig. U][b). The end charges and a 
background then disappear, so that the constraint of neu- 
trality becomes explicit, with conserving pi in form 
Next we decompose each charge above into fractions Ai , 
A 2 , and A3, which are then inverted relative to the ends 
of the unit cell in an antisymmetric fashion shown in Fig. 
Die). Taking into account the relation A\ + A 2 + A3 = 1 
normalizing the result per actual charges, the dipolar mo- 
ment pi of modified cell '1' can be written as 



Pi = 92 



12 



/21 



d k B! 



Pi + q2dkBi, 



(8) 



where B\ = A 2 — A3 is a new parameter. Thus the effect 
of inversion is realized in form (SJ as a shift of the q 2 - 
sublattice relative to the gi-one. 

The value of p 2 is defined in the same way, but with 
B 2 = —B± that supports the above treatment in terms 
of shifts and follows from the transformation p\ <-» p 2 
associated with (qi;h 2 ) <-> (92; hi)- Next we assume a 
functional form Bj — Dj/Dq so as to satisfy (0, where 
Dj are polynomials bilinear in qj and l^y , provided that 
Di = —D 2 is antisymmetric upon {q\\l\ 2 ) <-> (92^21), 
whereas Dq is symmetric upon both qi <-> q 2 and l\ 2 <-> 



3 



Z21 - Then the quantity of interest is of the form 

V,{ qi ^-,l,l) = ^ = - tX ^ lM , (9) 
D = d 2 k {ql +ql + hqiqs) + 2tl l2 l 2X qxq 2 , (10) 

with free parameters t and t\. Likewise, the same result 
arises upon considering cells '3' and '4', though i? 3 = 
— £>4 differs from B\ therein. It implies that the solution 
© (imp) is a mutual property of two sublattices, but not 
of their particular unit-cell representations. 

Here we tacitly imply that both q\ and q 2 are of the 
same sign and so homogeneous. In this case form (JZJ 
is also antisymmetric with respect to the replacement: 
qi — * h 2 , h 2 — * <72, <72 — > ?2i) and l 21 —* qi- In the event 
of @ and (fTUf> this fact results in t\ = 2. 

If (qiq 2 ) < 0, a similar consideration, but with \qiq 2 \ 
instead of qiq 2 in Dj, gives rise to ©, where 

D = d 2 k {ql +q\+ t 2 \ qi q 2 \) + 2tl 12 l 21 \ qi q 2 \. (11) 

In this case there is an additional limiting condition 



V 2 {q,-q;a,b} 



-aq 



(12) 



corresponding to a Bravais lattice of point dipoles with 
a period b and leading to t 2 = 2 — t. The parameter t is 
common to (|1U|) and from the continuity of @. On 
comparing the limit of @ as q\ — > oo and l\ 2 — > oo to 
JSJl, we obtain f = 2, because q 2 and Z21 then describe a 
small perturbation. Note that the other three choices of 
limiting couples of parameters are also admissible. As a 
result, we arrive at the normalized value 



V 2 {qi, q 2 ;h 2 ,hi) 

-l(qi,q 2 ;h 2 ,l 2 i) 



at {qiq 2 ) > 0, 



_ ) d 2 k ( qi + q 2 ) 2 + 4l 12 l 2iqi q 2 

-l(qi,q 2 ;h 2 ,l 2 i) , . 1 j 

rf|MTlP4W^ * (**,)< 0. 

The treatment above can be generalized to 71 sublat- 
tices described by an ordered set {qj} = {qi, . . . , q n } and 
by a complementary set {ljj+i} = {Z12, l 2 3, ■ ■ ■ , Ini} of 
nearest-neighbor distances. Keeping in mind that every 
couple of charges gives the independent contribution de- 
scribed by (|13fl . we readily reach 

1 " 

fe+i}) = - E E '^•h-'h - ! j. ■■!/,)■ (14) 

where + Zj'j = dfc for every j and j'. Every particular 
V 2 (qj,qj';ljj' ,lj' j) is counted twice in (JT2J, so that its 
own normalization is neutralized there. 

It is important that V n ({qj}; {ljj+i}) is not additive 
to the formal splitting of any nonzero qj into fractions, of 
which number would multiply the contribution of pj to 
V n {{qj}; {ljj+i}), with destroying the symmetry effect. 



Thus, every set of coinciding projections of point lattices 
must be reckoned once. 

To describe a real structure composed of m sublattices 
of point nuclei qj uc with a compensating electron density 
in the above manner, we map q^ uc onto the sets {qk,j} of 
nfc point charges, with retrieving the label k associated 
with efc, whereas the three electron charge projections 
pf(z) on these A-directions are defined like p(z) in 
Since any of pf(z) is consistent with the structural pa- 
rameters connected with {qkj} in equilibrium, 
its effect on (11411 is reduced to the replacement of qkj by 
effective charges: 



qk, 



Q 



i=i 



(15) 



where Q, the sum of qfj, is determined by the last equal- 
ity allowing for the total electrical neutrality. 

To determine qf . , a new set of unit cells based on the 

set {/j^} is to be introduced so as to avoid the degen- 
eracy associated with symmetric charge positions used 

(k) 

above. To this end, we note that any of the intervals 
is topologically equivalent, notwithstanding is it large or 
small relative to the others, and so must contain the 
boundary points of just one unit-cell species j of length 
dk- By homogeneity, the boundaries of the jth unit cells 
are assumed to be shifted by Ckl^fLi from the centers of 



the intervals ljj+i at a given j, with —1/2 < < 1/2 
common to all j. The values of q^. • are then determined 
by the set of rik linear equations which, at different j, 
describe the dipolar moments Fj k \^k) of the electron 
charge distributions within the jth unit cells in terms of 
the contributions of point charges q\ a, located properly: 



EC- 



Fl k) (£k 



(16) 



if to 



d k /2 
-d k /2 



z' P f\z k ,z')dz' = QRf\t k ). (17) 



Here the center of the jth unit cell is the origin common 
to the position r - k ) of the j'th charge, to the variable z' 



specifying pf\Ck,z') from pf(z), and to Rf'(^k), the 
center of gravity of the electron charge in that unit cell. 
On substituting the right equality of (|17fl into <|16|) . with 
taking 1)15(1 into account, the equations ((16JI become ho- 
mogeneous, so that their nontrivial solution exists at 



(fc), 



det 



33' 



= 0. 



(18) 



This is quite natural, because Q furnishes one more con- 
straint on qfj, . On casting rjy in terms of , one car 
prove by induction that relation l|18|) is converted into 



(fe) 



"A- 

7=1 



(fc) 

ii+i 



0. 



(19) 
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TABLE I: Pyroelectric vector and concomitant parameters of 
some ABO3 perovskites which are ferroelectric in a tetrago- 
nal phase, with model uniform spherical electron distributions 
employed. Experimental and Berry-phase-approach (BPA) 
results are also given for comparison. 





BaTi0 3 


PbTi0 3 


KNb0 3 


a (nm)° 


0.39947 


0.3904 


0.3997 


c (irni)" 


0.40336 


0.4150 


0.4063 


m a 


0.012 


0.041 


0.019 


m a 


-0.023 


0.112 


-0.025 


m a 


-0.014 


0.112 


-0.021 


t 


-0.0099 


0.0558 


-0.0329 


qx/e 


26.843 


31.093 


4.961 


qs/e 


0.944 


0.385 


14.774 


qi/e 


-22.416 


-24.280 


-5.434 


qil/e 


-5.371 


-7.198 


-14.301 


P° (Cm" 2 ) 


0.306 


-0.727 


0.358 


|Pe°x P l(Cm- 2 ) 6 


0.26 


0.75 


0.37 


I^PAKCm" 2 ) 


0.30; 0.28 c 


0.74; 1.04 d 


0.35; 0.40 e 




a BaTi0 3 and PbTi0 3 in Ref. E3, KNb0 3 in Ref. |24|. 
^BaTiC^ in Ref. E|, PbTi0 3 in Ref. KNb0 3 in Ref. 
and Ref. [g, respectively. 

calculations at 700 K and at 295 K, respectively! 
and Ref. 0, respectively. 



With making use of the left equality of l|17(l . equation 
(TH^I specifies a self-consistent value of while set (fLo)) 
gives rise to qfj, , so that q~k,j ln nas been definite. 

On joining all the results obtained together, the pyro- 
electric density vector is well-defined by the form 



p0 = ^E e ^({^}'{4 fc li})' 



(20) 



fc=l 



where v is the unit-cell volume. Deduced from pristine 
charges by linear operations, l|2U|) agrees with 

Note that in the phase-transition problem relation l|20() 
implies that a sole nonzero component of P° can arise 
without jumps, with treating P° as an order parameter 
small near the transition point Furthermore, 
different nonzero projections of P° in (|2U[) can arise sep- 
arately, as it happens in BaTiC>3 and KNbC>3 [2! EH |22| . 

Concentrating on regular states, here we estimate the 
polarization of BaTiC>3, PbTi03, and KM3O3 in a tetrag- 
onal phase. In general terms of ABO3 perovskite, such a 
ferroelectric phase is determined by the lattice constants 
a and c and by extra displacements ctjb, crji, and cijn of 
the B, Oi, and two On sublattices relative to the A one, 
provided that Oi and On ions surround B along c and 
in lateral directions, respectively, and positive 7/b defines 
the polar direction. The simple model of uniform electron 
distributions within ionic spheres of Ba 2+ , Pb 2+ , Ti 4+ , 
K + , Nb 5+ , and O 2- ions, with On anions in contact, is 
employed. With making use of the structural parameters 



listed in Tableland with account of the atomic numbers 
Z Ba = 56, Z Ti = 22, Z = 8, Z Ph = 82, Z K = 19, 
and -ZNb = 41, the effective charges qj are obtained as 
a solution of equations H15|l -119 (l and are listed in Table 
[I] in units of the elementary charge e. Then the values 
of P° are obtained with the help of JT3J), ifHjl. and (|2U|) 
and agree with the experimental data, as shown in Ta- 
ble |U On comparing with the results of the Berry-phase 
approach shown in Table [fl as well, our estimate is much 
closer to the experimental data on PbTi03 at room tem- 
perature |26j |. but is comparable in other cases. 
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